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Abstract 

We consider the low-energy electronic properties of graphene cones in the presence of a global 
Fries-Kekule Peierls distortion. Such cones occur in fullerenes as the geometric response to the 
disclination associated with pentagon rings. It is well known that the long-range effect of the 
disclination deficit-angle can be modelled in the continuum Dirac-equation approximation by a 
spin connection and a non-abelian gauge field. We show here that to understand the bound states 
localized in the vicinity of a pair of pentagons one must, in addition to the long-range topological 
effects of the curvature and gauge flux, consider the effect the short-range lattice disruption near 
the defect. In particular, the radial Dirac equation for the lowest angular-momentum channel sees 
the defect as a singular endpoint at the origin, and the resulting operator possesses deficiency 
indices (2,2). The radial equation therefore admits a four-parameter set of self-adjoint boundary 
conditions. The values of the four parameters depend on how the pentagons are distributed and 
determine whether or not there are zero modes or other bound states. 

PACS numbers: 73.63.-b; 05.30.Pr; 05.50.+q 
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I. INTRODUCTION 



The essential features of the valance and conduction bands of planar graphene [1] can 
be modelled as a pair of 2+1 dimensional Dirac fermions, one for each of the conical band 
touchings that occur at the points k = K and k = K' in the Brillouin zone [2]. A Fries [3] 
(or equivalently Clar [4]) Kekule-structure distortion of the hexagonal graphene lattice will 
couple the two Dirac points and introduce a mass gap. If the phase of the K-K' coupling can 
vary with position, there may exist topologically stable vortex textures that are associated 
with zero modes and charge fractionalization [5-7]. 

A complex-valued K-K' coupling is unlikely to occur spontaneously as a result of a sim- 
ple Peierls distortion - although it might be artificially engineered through proximity-effect 
coupling to vortices in a superconducting substrate [7, 8]. An alternative and naturally oc- 
curring form of vortex may be provided by local curvature-inducing geometric defects such 
as pentagons and heptagons [9-12]. A spherical fullerene such as C 60 possesses twelve pen- 
tagonal defects and, even in the absence of Kekule distortion, its energy spectrum contains 
six delocalized low-energy levels. These near-zero-energy levels have long been understood 
as the lattice relicts of six exactly zero-energy modes of the continuum Dirac hamiltonian 
[13, 14]. The continuum zero-modes are those predicted by the index theorem for 2+1 di- 
mensional Dirac fermions moving in a fictitious monopole magnetic field that mimics the 
effect of the defects. One quarter of a Dirac unit of monopole flux threads through each 
of the twelve pentagons, and for Ceo this discretely-lumped field is sufficiently spread out 
that the resulting energy spectrum is well approximated by a spatially uniform field [13, 14]. 
There are three units of flux altogether, and so the index theorem predicts three zero modes 
apiece for the two Dirac fermion species of planar graphene. 

The fictitious flux through the closed surface of the fullerene molecule means that when 
the discrete lattice hamiltonian is approximated by two continuous-space Dirac hamiltonians, 
the continuum wavefunctions must be regarded as sections of a twisted line bundle, the 
twist being characterized by a Chern number of ±3. When a purely real Kekule field is 
introduced on the fullerene, the threefold twist, when coupled with a natural choice of 
how we introduce the gauge holonomy, allows the real K-K' coupling matrix elements to 
be perceived by the continuum approximation as being associated with a charge-2 complex 
Higgs field possessing a net winding number of six. Each pair of pentagons is then effectively 
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a single vortex, and each such vortex should, by the Jackiw-Rossi- Weinberg index theorem 
[15, 16], contain at least one zero mode. The six low-lying "zero modes" should therefore 
survive the introduction of the Kekule induced mass gap, and, in theory, its principal effect 
should be to localize the previously extended states in the vicinity of the pentagonal defects. 

In practice, numerical investigation of the lattice spectrum [11] shows that although the 
"zero modes" are not immediately destroyed by the introduction of the Kekule Higgs field, 
their energy is affected. Indeed as the ratio h of the double bond to single bond hopping 
evolves from zero to large values the "zero-modes" cross the finite-size gap from the negative 
to the positive part of the spectrum. What was not obvious from the plots in [11] is that 
at the same time the zero-mode wavefunctions evolve from being tightly localized around 
the defects to being antilocalized . Thus these modes, while still topologically interesting, 
are not acting as expected from the continuum Jackiw-Rossi- Weinberg index theorem. This 
theorem predicts that the energy of the zero modes will be unaltered by the introduction of 
the Higgs field, and that they will be localized for any sign of the mass. 

In this paper we will explore what feature of the continuum approximation to the lattice 
hamiltonian is responsible for this behaviour. We will focus on isolated pentagons and 
isolated pairs of pentagons. These disinclination defects roll the sheet of graphene into a 
cone. Away from the tip of the cone the geometry is locally flat and a continuum twisted 
Dirac hamiltonian approximation should be reliable. Near the tip, the exact form of the 
lattice becomes important. The lattice effects can be accounted for, however, by imposing 
a boundary condition on the continuum wavefunction at the tip of the cone. The result is 
the introduction of a four-parameter family of self-adjoint extensions to the Dirac operator. 
General values of the extension parameters explicitly break the symmetry required for the 
index theorem, and suitable choices of the parameters reproduce the numerically observed 
phenomena. 

II. PLANAR GRAPHENE AND THE DIRAC EQUATION 

In order to establish our notation, we begin with a quick review of how the tight-binding 
(Hiickel) approximation [13] for planar graphene leads to a massive Dirac hamiltonian. In 
this approximation the 7r electrons hop on a two-dimensional honeycomb lattice with lattice 
constant a. The honeycomb lattice is bipartite with the property that nearest-neighbour 
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FIG. 1: The zero-energy reference states \A+), \B+), \A—) and \B— ). The empty circles indicate 
that the wavefunction is zero at those sites. The symbols rj and tj 2 indicate that the wavefunction 
takes the values 77 = exp(27ri/3) and if = fj = exp(— 27ri/3) at that site. The heavy lines indicate 
a Fries stucture of "double" bonds, whose hopping will increase from t to t + St when a Peierls 
distortion occurs. 

hopping takes an electron from sublattice A to sublattice B, and vice versa. We denote the 
hopping amplitude between nearest neighbour sites in undeformed graphene by t. A Kekule 
structure is an assignment to the lattice edges of a pattern of single and double bonds 
such that each carbon atom partakes of two single bonds and one double bond. There are in 
general many such assignments, but patterns that lead to a maximum number of benzene-like 
hexagons are known as Fries structures. Fullerenes that possess a globally defect-free Fries 
structure (the so-called leapfrog fullerenes [17]) are typically the most chemically stable. 
This is because a spontaneous Peierls distortion that shortens the Fries-structure double 
bonds and increases their hopping amplitude from t to t + 8t will open a gap between the 
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highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital 
(LUMO), lower the energy of the system, and give the molecule the property of having a 
"closed shell." Such a spontaneous distortion occurs in C 60 where the double bonds have 
length ~ 0.1388 nm while the single bonds have length ~ 0.1432 nm [18]. 

The hopping hamiltonian for the undistorted infinite lattice has four linearly independent 
zero-energy eigenstates which are displayed in figure 1. The labels A and B indicate that 
the non-zero amplitudes are supported on sublattice A or B respectively, and the "+" and 
'— " configurations have opposite handedness in the way that their complex phases evolve as 
we circle the zero sites. Both "+" configurations have the same lattice momentum k = K 
and both "— " configurations have lattice momentum K'. Any low energy eigenfunction will 
be a slowly varying linear combination 
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of these reference wavefunctions. Here r labels discrete lattice points, but since the /(r) 
are to be slowly varying on the scale of the lattice they can be regarded as being smooth 
continuum functions M 2 — > C. 

The hamiltonian acts on the array of slowly varying functions / as 
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When the hopping is modified by the Fries-Kekule structure shown in figure 1, the reference 
states are no longer exactly zero-energy eigenfunctions. Instead H\A+) = St\B—) and 
H\A—) = 5t\B+). Non-zero matrix elements therefore appear coupling the K Dirac point 
to the K' point: 
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where $ = (2/3at)5t. The derivative blocks can be written as 
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where cr = (<t 2 , — Ci) is a 90° rotated version of the cr = (a±, <7 2 ) spin vector. 
We can if, we prefer, write 
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in which case the derivative blocks involve the more familiar form 
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We thus reduce the hamiltonian to an unimportant constant times 
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This is a two-dimensional reduction of the mass m = |$|, three-dimensional Dirac hamilto- 
nian. We notice that the matrix 

V 3 o 

-<T 3 

anti-commutes with ifoirac- This reflects the fact that acting on a wavefunction by T inverts 
the sign of the wavefunction on one of the two sublattices whilst leaving the other sublattice 
alone, and that this involution takes an eigenfunction of energy E to one of — E. Zero energy 
states can be chosen so that they are eigenvectors of T, and the number of +1 eigenvalues 
minus the number of —1 eigenvalues is related to topological data by the Jackiw-Rossi- 
Weinberg index theorem. 



III. GEOMETRIC DEFECTS AND THEIR GAUGE FIELD 

We insert a pentagon disclination into the lattice by cutting out a 60° wedge and rejoining 
the severed bonds. We wish to preserve the global Fries-Kekule pattern, so the apex of the 
removed wedge must lie in a hexagon with no double bonds (see figure 2). In order for the 
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FIG. 2: Cutting out a 60° wedge and reconnecting the severed bonds leaves a pentagon disclination. 
When the apex of the wedge lies in a hexagon with no double bonds, single bonds reconnect to single 
bonds and double bonds to double bonds, leaving a globally consistent Fries-Kekule structure. 



lattice wavefunction to join smoothly across the seam of new bonds, the coefficient functions 
on the original lattice must obey the boundary conditions [9, 10] 

ff A+ \ /0 fj 0\ ff A+ \ 
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Here the subscript R means the value to the right of the cut, and L to the left. In order to deal 
with this transformation it is convenient to change basis so that \& = (/a+> /b+, /b-, — /a-) t - 
This makes ± kinetic energy blocks become identical 
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at the expense of complicating the mass terms. We can write the Hamiltonian matrix 
compactly as 

H= (I® or) •p + $ri®cr 3 , (12) 



where the I and T\ matrices act on the two-by-two blocks (i.e. in the K-K' space) and the 
a matrices act between the A and B sublattices within each K, K' block . (In the sequel 
we will omit I matrices when no ambiguity results.) 

In the new basis, the operator V that anti-commutes with H remains 

r = r 3 <g> 0-3 = 

The boundary condition, however, becomes 
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This can be factored as 



\&r = — iT\ <S> exp{— 

= exp{— inTi/2} ® exp{— ma^/Q}^^. 



(14) 



As before, the T\ matrix acts on the ± valley-degeneracy "flavour" indices and the a matrices 
on the A, B "spin" indices. 

Similarly, cutting out a 120° wedge turns a hexagon into a square, and requires 
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which can be factored as 



^ R = -I® exp{-ina 3 /3}^ L 

= exp{—inTi}®exp{—i'Ko- 3 /3}^L. 



(16) 



In this last line, the T\ matrix in the exponent can be replaced by r 2 , r 3 , or even by I, without 
altering the lattice boundary condition. The authors of [11] elected to take this matrix to 
be T3 as this choice leads to a continuum hamiltonian for which the symmetry required by 
the Jackiw- Rossi- Weinberg index theorem appears manifest. 
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We can remove the discontinuity across the 60° wedge by writing 

tt(r,0) =exp li^-n ( |-U ®exp \£a 3 (f^ j W(r,0), (17) 
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and across the 120° wedge by writing 

tt(r, 0) = exp jmn J ® exp ji|<7 3 | *(r, 0). (18) 

Again, in the second case, we may replace the T\ matrix by r 2 , r 3 or I. In all cases, the 
new field ^(r, 0) is continuous across the reconnected seam. In the first case the angle is 
restricted to the range — n/6 < 9 < 3ix/2, with the limiting values representing the same 
point on the graphene cone, and in the second n/6 < 9 < 3ii/2. 

We wish to write the eigenvalue problem in polar coordinates whose origin is at the tip 
of the cone. To do this we use the identity 
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The coefficients of the derivatives a r> e = 61,2 are now matrix-valued constants, but they 
pay for their constancy by being attached to a moving zweibein frame e r , e#. The — icr 3 /2 
spin connection is therefore required to cancel the effect of taking the frame-rotation matrix 
exp {ia 3 0/2} through the 9 derivative. 

For planar graphene we would now define a new field x( r ; 9) = e^ a3d ^(r, 9) and find that 
it is antiperiodic under 0^0 + 2n. Because of the deleted wedge, however, we have to 
define 

X (r, 9) = exp j«7 S Q + 1) J *(r, 0), (20) 

which is antiperiodic under 0^0 + 27r(5/6). Taking into account eq. (17), the eigenvalue 
problem H^> = E^fr for the 60° wedge can now be written 
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Here the «7~i(3/10) gauge connection comes from taking the exp {ir\ (3/10)0} matrix appear- 
ing in equation (17) through the derivative. For the 120° wedge we must set, 
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which is antiperiodic under 9 — > 9-\-2ti (2/3). The corresponding eignevalue problem becomes 
.„ d . 1 
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From now on, for purely cosmetic reasons, we will make the rotation introduced in (6) 
that removes the tildes from the a\^i matrices. This does not affect the gauge fields, as they 
act in the valley degeneracy "flavour" space and the redefinition acts in the A-B "spin" 
space. 

We define a new angle = (6/5)9, or = (3/2)# that has the usual 2ir periodicity. (0 is 
the polar angle seen when looking along the axis of the 60° or 120° cones from above their 
apex.) Using we can then separate the radial and angular part of the wavefunction as 
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Note that the hermiticity of the radial part of the differential operator with respect to 
the inner product 

^ 2 rdrd(t) (26) 

requires a contribution from the spin-connection term —203/2 that occurs in the angular 
covariant derivative part: 
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Because T\ commutes with H, we can find solutions to (21) and (23) (without the tilde's) 
of the form 
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which are eigenvectors of ri with eigenvalue ±1. The functions a(r),v(r) then satisfy 
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Here n = 1 for the 60° cone and n = 2 for the 120° cone. The number r = ±1 denotes the 
eigenvalue of t±. Set 
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We have continuous-spectrum eigenfunctions with u, v of the form 
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The first set of solutions is finite at the origin when j > 1/2 and the second is finite at 
the origin when j < —1/2. For j — 1/2 and r negative, the upper component of the first 
set of solutions diverges at the origin, but no faster than r _1//2 , so it is normalizable there. 
Similarly the second set is locally normalizable for j = —1/2. 

For most values of j and r, demanding nomalizability is sufficient to select the physically 
allowed solutions. When n — 2, however, and for j = 1/2, r = —1, we have v — 0. We also 
have v — when j = —1/2, r = +1. In these cases, the eigenvalue equation becomes 
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The scattering solutions with _E = ±V^ 2 + $ 2 contain the Bessel function Ji/ 2 {kr) = 
sin/cr and J-i/ 2 {kr) = ^2n/kr cos A;r, both of which are normalizable near the 
singular point at the origin. The equation is therefore in Weyl's limit-circle class there, 
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and some additional boundary condition must be imposed to select a complete, linearly 
independent, set of solutions [19]. 

One reason for the n — 2, 120° wedge differing from the n — 1, 60° wedge is that that 
the cutting and sewing of the graphene sheet in the former case preserves the A-B bipartite 
structure of the lattice. There should therefore be a some operator that anti-commutes with 
the Hamiltonian. This property is not easy to see in eq (23), but becomes clearer if we elect 
to replace the gauge field term with r 3 . Then equation (18) defining the single-valued field 
^(r, 9) becomes 
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with antiperiodic x( r >@)- This is the equation with a unit- winding-number vortex in the 
mass term that was considered in [11]. It is formally of the form to which the Jackiw- 
Rossi- Weinberg index theorem applies: the matrix-valued differential operator manifestly 
anti-commutes with r = r 3 <g> <r 3 , and the index is the Hilbert-space trace of T. 
Eq (36) posseses zero energy solutions of the form 
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which are eigenvectors of T with eigenvalues +1 and —1 respectively. Ignoring any con- 
straints imposed by boundary conditions they are 
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Only one pair will be normalizable, depending on the sign of <£>. 

Eq (36) also posses more general bound-state solutions for any value of E in the range 
-|$| < E < |$|. These are 
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Here /t = \f§ 2 — -E 2 . These solutions are not eigenf unctions of T. Instead T acts on them 
to give a solution with energy — E. 

Both the E = and the more general — 1$| < E < |$| solutions are square integrable 
in the neighbourhood of r = 0. The same is true of the scattering state solutions with 
\E\ > |$|. We need to impose some boundary condition at r = to select from these 
solutions a complete linearly-independent set. As a guide to the form of this boundary con- 
dition we apply the Weyl-von-Neumann theory [20] . The matrix- valued differential operator 
determining the radial part of all these solutions is 



H = —io\ 



(42) 



If we restrict it an initial domain that contains only functions that vanish at r = 0, the 
resulting linear operator has deficiency indices (2,2). The restricted operator therefore 
admits a four-parameter family of self-adjoint extensions. To determine the corresponding 
boundary conditions, we evaluate (&\HX) — (H^\X) on functions ^ = (ipi, ip 2 , ^3, i>A) T and 
X = (xi, X2, X3? Xi) T that are square integrable on [0, oo) with the measure rdr. We find 
that 

($\HX) - {H*\X) = [-1^1X2 + ^2X1 + + r,X3)\o- (43) 

Since \J/(r) and X(r) are allowed to diverge as r" 1 / 2 near the origin but must tend to 
zero faster than r~ 2 at infinity, the vanishing of the integrated-out part requires that the 
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FIG. 3: The left-hand plot shows F(E) for boundary-condition parameters A = B = D = 0, C = 1, 
and $ = — 1 which possesses two degenerate zero modes. The right-hand plot shows F(E) for the 
same boundary parameters, but with 3> having changed sign from negative (reduced hopping in 
the double bonds) to positive (enhanced hopping on the double bonds). There are now no bound 
states. 







expression in parentheses tend to zero at r = 0. If we impose boundary conditions 
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on \l/(r), the adjoint boundary conditions on X{r) are determined by requiring that 
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for any ■02? i ) 4- Thus the adjoint boundary conditions are 
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For H to be self-adjoint, these boundary conditions must coincide with the boundary con- 
ditions imposed on ip. We there have that, as r — > 0, 

VhA / iA B + iC\ U 2 
^) ~ \-B + iC iD ) \> 

for real (possiby infinite) numbers A, B, C and D. This relation contains the four real 
parameters required by the Weyl-von Neumann count, and so is the most general possible 
self-adjoint boundary condition. 

The numerical solution of the 120° wedge-cut lattice show that there are two exact zero 
modes that are localized when $ < and become delocalized when $ > 0. These modes 



14 



A=0.0, B=0.<) 
C=1.0, D=0.5 
=-1.0 



REi 



A=-0.S, B=0.0 
C=l.l,D=-0.5 
0.8 - «=-1.0 




FIG. 4: The left-hand figure shows F(E) for A = B = 0, C = 1, D = .5 and $ negative. There is a 
zero mode and a bound state with E ~ — 0.5|<3?|. In the right-hand figure we have A = D = —1.0, 
B = 0, D = y/2, and negative. There are two degenerate bound states at E = |<fr|/\/2. In both 
cases, the bound states cease to exist as soon as <!> changes sign from negative to positive. 



(48) 



can be identified with ^o++ and ^o— h an d are allowed if we impose the boundary condition 

>i\ _ /o A 

at r = 0. The other potential zero modes $oh and do not satisfy this condition. 

For general values of A, B, C and D, we can seek bound-state solutions of the form 

Imposing the boundary condition (47) at r = leads to a pair of homogeneous equations 



{a + p)E + {a-/3)$ 
{a-{3)E + (a + [3)$ 



-An{a + /3)-(C- iB) K {a - p), 
-Dn{a -P)-(C + iB)n{a + p), 



(49) 



for a and P, and hence to the condition F(E) = 0, where 

E + Ak <&+(C-iB)K 

F{E) 



(50) 



§ + (C + iB)K E + Dk 

The function F(E) is real in the range — 1$| < E < |$|, and always has two zeros at E — 
±|$| corresponding to the edges of the upper and lower continuum respectively. Additional 
zeros in the range — 1$| < E < |$| correspond to the energies of bound states. Example 
plots of F(E) for different values of A, C, D and $ are shown in figures 3 and 4. Of 
particular interest is the left hand plot in figure 3 which exhibits the pair of E = zero 
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FIG. 5: The low- lying part of the energy spectrum plotted versus h (the ratio of double bond 
hopping to single bond hopping) for a square defect created by excising a 120° wedge. The hori- 
zontal axis has been displaced vertically so as to uncover the doubly-degenerate exact zero mode. 
The "within-gap" modes peeling off from the upper and lower continua for h > 1 {i.e. $ > 0) are 
"anti-localized" edge states, whose exact form depends on how we truncate the lattice on its outer 
boundary. 

modes for the boundary conditions in (48), and the right-hand plot in figure 4 which exhibits 
a pair of degenerate levels at a non-zero value of E. The plots in figure 4 correspond to 
non-zero values for one or both of the parameters A and D. These parameters control the 
boundary coupling of the A and B lattices to themselves, and non-zero values violate the 
bipartite lattice structure that leads to the E <-» — E spectral symmetry possessed by the 
bulk differential equation. 

Figures 5, 6 and 7 show some numerical plots of the low- lying energy states as a function 
of h = (t + 5t) 1 1, where St is the change in hopping parameter on the double bonds (recall 
that $ = (2/3at)5t.). In all four lattices, the wedges have their apices in single-bond 
hexagons, so that the global Kekule structure is preserved. In the language of [9] they all 
have n = mmod3, ensuring that the asymptotic gauge field is the same in all four cases. 
The spectra differ, however, because the boundary conditions seen by the continuum wave 
functions at the tip of the cone are different. In the case of figure 5, the global bipartite 
structure is preserved, the boundary conditions are those of eq. (48) and so the spectrum 
is manifestly E <-> —E symmetric. In lattices of figures 6 and 7 the bipartite structure is 
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FIG. 6: The low-lying energy spectrum plotted versus h for a pair of nearby pentagons. The 
left-hand figure is for an (n.m) = (1,1) cone in the languge of ref [9], and the right-hand figure is 
for an (n, m) = (0, 3) cone. In both cases the 60° wedges have their apices in single-bond hexagons 
so as to preserve the global Kekule structure. For h < 1 there is a pair of nearly degenerate bound 
states lying just below the upper continuum. The "below gap" modes at h > 1 are uninteresting 
edge states localized at the outer boundary. 

scrambled in the region between the pentagons and so the E <-> —E symmetry is violated. 
All of the two-pentagon cases display a pair of almost degenerate bound states just below the 
positive continuum, and are therefore similar to the spectrum associated with the boundary 
conditions of the right-hand plot of figure 4. As h approaches unity, and $ approaches zero, 
the wave functions begin to spread out and see the outer boundary of our large-but-finite 
lattices (2644 vertices in case of the square and 1158 vertices for the pair of pentagons). 
The effects of this can be seen in the figures in the splitting of the bound-state energies as 
h approaches unity from below. 

IV. CONCLUSIONS 

The continuum Dirac hamiltonian provides a good account of the low-energy, long- 
wavelength, eigenstates of the tight-binding hamiltonian on an infinite sheet of graphene. 
The continuum model is also a useful approximation for cones tipped by curvature singular- 
ities induced by pentagon defects — but it must be supplemented by non-trivial boundary 
conditions at the tip of the cone. Although the low-lying eigenfunctions have too long a 
wavelength to resolve the fine details of lattice disruption at the tip, they there experi- 
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FIG. 7: The low-lying energy spectrum plotted versus h for a pair of pentagons forming a (n, m) = 
(0, 6) cone. 

ence phase shifts and mode mixing that have a significant effect on the eigenstates. For 
two separated pentagons, the scrambling of the A-B bipartite lattice structure along a seam 
joining the pentagons sufficiently violates the E <-> —E spectral symmetry as to allow bound 
states at non-zero E. In the case of two coincident pentagons {i.e. a square) the E <-> —E 
symmetry is preserved, but (in contrast to the quarter-unit flux through a pentagon) the 
resultant half-unit of flux through the square plaquette is too large to be approximated by 
a spread-out gauge field. A continuum gauge field with this flux would have bound a single 
state whose eigenvalue of Y is determined by the sign of the flux. The lattice spectrum must 
be unchanged, however, by the insertion of an integer flux-quantum through the square. 
Such an insertion can reverse the sign of the flux and so no particular sign of Y can be 
favoured. The lattice hamiltonian compromises by producing two bound states, one with 
each sign of Y . All these effects can be reproduced in the continuum model by suitable 
choices of the parameters in the self-adjoint boundary conditions. However, the predictions 
[11] of the continuum Jackiw- Rossi- Weinberg theorem do not survive in a simple form. It 
will be interesting to see if the theorem can be modified to include the effects of the singular 
endpoint. 
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